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ABSTRACT 

We investigate particle accelerators in rotating neutron-star magnetospheres, by simultane- 
ously solving the Poisson equation for the electrostatic potential together with the Boltzmann 
equations for electrons, positrons and photons on the poloidal plane. Applying the scheme to the 
three pulsars, Crab, Vela and PSR B1951+32, we demonstrate that the observed phase-averaged 
spectra are basically reproduced from infrared to very high energies. It is found that the Vela's 
spectrum in 10-50 GeV is sensitive to the three-dimensional magnetic field configuration near the 
light cylinder; thus, a careful argument is required to discriminate the inner-gap and outer-gap 
emissions using a gamma-ray telescope like GLAST. It is also found that PSR B1951-I-32 has a 
large inverse-Compton flux in TeV energies, which is to be detected by ground-based air Cerenkov 
telescopes as a pulsed emission. 

Subject tieadings: gamma-rays: theory - magnetic fields - methods: numerical - pulsars: individ- 
ual(B1951+32,Crab,Vela) 



1. Introduction 

Pulsars form the second most numerous class of 
objects detected in high-energy 7-rays. To date, 
six have been detected by the Energetic Gamma 
Ray Experiment Telescope (EGRET) aboard the 
Compton Gamma Ray Observatory. The 7-ray 
pulsations observed from these objects are par- 
ticularly important as a direct signature of non- 
thermal processes in rotating neutron-star magne- 
tospheres, and potentially should help to discrim- 
inate among different emission models. 

The pulsar magnetosphere can be divided into 
two zones: The closed zone filled with a dense 
plasma corotating with the star, and the open zone 
in which plasma flows along the open field lines to 
escape through the light cylinder. The last-open 
field lines form the border of the open magnetic 
field line bundle. In all the pulsar emission mod- 
els, particle acceleration takes place in this open 
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zone. In inner-gap (IG) models, which adopts par- 
ticle acceleration within several neutron-star radii 
above the polar-cap surface, the energetics and 
pair cascade spectrum have had success in repro- 
ducing the observations (e.g., Daugherty & Hard- 
ing 1982, 1996). However, the predicted beam 
size is too small to produce the wide pulse pro- 
files that are observed from high-energy pulsars. 
Seeking the possibility of a wide hollow cone emis- 
sion due to flaring field lines, Muslimov and Hard- 
ing (2004) extended the idea of Arons (1983) and 
proposed a slot-gap (SG) model, in which emis- 
sion takes place very close to the last-open field 
line from the stellar surface to the outer magneto- 
sphere. Since the SG model is an extension of the 
IG model into the outer magnetosphere, a nega- 
tive magnetic-field-aligned electric field, , arises 
if the magnetic moment vector points in the same 
hemisphere as the rotation vector. However, the 
electric current induced by the negative iSy does 
not have a self-consistent closure within the model 
(Hirotani 2006, hereafter H06). 

To contrive an accelerator model that predicts 



an outward current in the lower latitudes (within 
the open zone), it is straightforward to extend 
outer-gap (OG) models (Cheng et al. 1986; Ro- 
mani & Yadigaroglu 1995, hereafter RY95; Cheng 
et al. 2000) into the inner magnetosphere. Ex- 
tending several OG models (Hirotani et al. 2003; 
Takata et al. 2004), HOG demonstrated that the 
gap extends from the stellar surface to the outer 
magnetosphere, that the positive extracts 
ions from the star as a space-charge-limited flow 
(SCLF), and that most photon emission takes 
place in the outer magnetosphere because E\\ is 
highly screened inside the null surface owing to the 
discharge of the created pairs. In the present let- 
ter, we formulate the scheme (Beskin et al. 1992) 
in § [21 apply it to three rotation-powered pulsars 
in § [3l and give a brief discussion in § [4l 

2. Gap Electrodynamics 

We follow the scheme described in § 2 of H06 to 
solve the set of Maxwell and Boltzmann equations. 
The first equation we have to consider is the Pois- 
son equation for the electro-static potential, 5*. As 
space is limited, we present its Newtonian expres- 
sion, — V^^* = 47r(p — pgj)- If the real charge den- 
sity, p, deviates from the Goldreich- Julian (GJ) 
charge density, pc], in some region, an acceler- 
ation electric field E'y = —d'^/ds arises, where s 
designates the distance along a magnetic field line. 

The second equation we have to consider is the 
Boltzmann equations for e^'s. Assuming a steady 
state in the frame of reference corotating with the 
magnetosphere, we obtain 
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where c denotes the speed of light, n+ (or n_) the 
dimensionless positronic (or electronic) distribu- 
tion function normalized by the local GJ number 
density. The temporal derivatives of momentum 
and pitch angle, dp/dt and dx/dt, and the colli- 
sion term S± are explicitly given in HOG. Synchro- 
curvature radiation-reaction force is included as an 
external force in dp/dt, while the effects of inverse- 
Compton scatterings (ICS) and (one-photon and 
two-photon) pair creation processes are in S±. 

The third equation we have to consider is the 
Boltzmann equations for photons. Assuming a 
steady state, and neglecting azimuthal propaga- 



tions, we obtain 

where the wave numbers fc*" and are given by 
the ray path, {r,0) are the polar coordinates, and 
the dimensionless photon distribution function g is 
normalized by the GJ number density at the stel- 
lar surface. ICS, synchro-curvature emission, and 
the absorption are contained 111 S'y . In H06, the 
rate of synchrotron emission by secondary pairs 
created outside the gap, was calculated assuming 
a constant pitch angle. However, it turns out that 
only 17 % of the initial particle energy is con- 
verted into radiation. In this letter, we corrected 
this problem by computing the pitch angle evo- 
lution of radiating particles, which increases the 
synchrotron cooling time and hence recovers the 
time-integrated, radiated energy to 100 % of the 
initial particle energy. Note that the gap electro- 
dynamics investigated in H06 remains correct, de- 
spite the insufficient secondary synchrotron fluxes 
in H06. 

To solve the Poisson equation, we impose = 
at the lower, upper, and inner (s = 0) boundaries, 
and — —d'^/ds = at the outer boundary. 
Ion extraction rate is regulated by the condition 
i^ll = at s = 0. We parameterize the trans-field 
thickness of the gap with h„^. If = 1 the gap 
exists along all the open field lines, while if ft-m ^ 1 
the gap becomes transversely thin. 

3. Application to Individual Pulsars 

We apply the theory to three 7-ray pulsars. 
Crab, Vela, and B1951-I-32, focusing on their pho- 
ton spectra. Even near and outside the light cylin- 
der, photon emission and absorption occur effec- 
tively; thus, equations H]) and ^ are solved in 
< s < 16n7Lc, where E\\ = in -cu > 0.9n7Lc; 
tULc = c/r2 denotes the light cylinder radius, f2 
the stellar rotation frequency, and vu the distance 
from the rotation axis. The field line geometry in 
0.9'n7LC < vj < 2ci7LC mimics the aligned dipole 
in the force-free magnetosphere (Contopoulos et 
al. 1999; Gruzinov 2005). In to > 2tnLc, the field 
lines are assumed to be straight so that they con- 
nect smoothly at zu = 2wi^c- 
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3.1. Crab pulsar 



3.2. Vela pulsar 



We present the results for the Crab pulsar in fig- 
ure [U adopting a magnetic inclination of ai = 75° 
and a dipole moment of ^ = 4 x 10'^° G cm'^, which 
is close to the value (3.8 x 10'^°Gcm^) deduced 
from the dipole radiation formula. It follows that 
the observed pulsed spectrum from IR to VHE 
can be reproduced, provided that we observe the 
photons emitted into 75° < 6 < 103°, where de- 
notes the photon propagation direction measured 
from the rotational axis. Because of the aberra- 
tion of light, it is reasonable to suppose that pho- 
tons emitted in a certain range of 9 comes into 
our line of sight in an obliquely rotating three- 
dimensional magnetosphere (e.g., RY95). The flux 
rapidly decreases with decreasing 9 for 75° < 9 < 
93°, because -Ey is highly screened in the inner 
part of the gap. Nevertheless, an average over 
75° < 9 < 103°, which includes neghgible flux 
between 75° < 9 < 89°, achieves the current ob- 
jective, because the spectral normalization can be 
fitted (within a factor of a few) without changing 
the spectral shape, by adjusting h„^. 

ICS takes place efficiently near and outside the 
light cylinder (Aharonian & Bogovalov 2003), be- 
cause the magnetospheric IR photons, which are 
emitted along convex field lines, collide with the 
gap-accelerated positrons near the light cylinder, 
where the field lines are concave. Most of such 
upscattered photons, as well as some of the high- 
energy tail of the curvature component, are ab- 
sorbed by the 7-7 collisions. As a result, there is 
a gradual turnover around 10 GeV, which forms 
a striking contrast with the steep turnover pre- 
dicted in IG models due to magnetic pair creation. 
The primary curvature component appears be- 
tween 100 MeV and 10 GeV, while the secondary 
synchrotron component appears below a few MeV. 
Between a few MeV and 100 MeV, the ICS compo- 
nent dominates, because the secondary pairs that 
have been cooled down below a few hundred MeV 
efficiently up-scatter magnetospheric UV and X- 
rays to lose energy. Similar spectral shapes are 
obtained for different values of cti, ft-m, provided 
that the created current is super GJ, by virtue of 
the negative feedback effects (II06). 



Next, we present the results for the Vela pul- 
sar in figure [5] (left). Taking an angle average 
over 75° < 9 < 103° (solid fine), we can repro- 
duce the observed pulsed spectrum, except for the 
RXTE results. The primary curvature component 
appears between 100 MeV and 10 GeV, while the 
secondary and tertiary synchrotron components 
appear below 100 MeV. ICS is negligible for the 
Vela pulsar because of its weak magnetospheric 
emission. Similar spectral shapes are obtained for 
super-GJ solutions, even though we have to adjust 
/im to obtain an appropriate flux. 

In the right panel, we compare the present re- 
sults with IG (dotted) and OG (dashed) models, 
where the dash-dotted line denotes the averaged 
flux for 75° < 6* < 107°. It follows that the spec- 
trum between 10 and 100 GeV crucially depends 
on the angles in which the photons that we ob- 
serve are emitted. Thus, to quantitatively predict 
the 7-ray emission from the outer magnetosphere, 
it is essential to examine the three-dimensional 
magnetic field structure near and outside the light 
cylinder. 

3.3. PSR B1951+32 

Thirdly and finally, we present the spectrum 
of B1951-h32 in figure [3] (left). It follows that 
the ROSAT and EGRET fluxes are reproduced by 
taking the flux average in 75° < 9 < 103° (solid 
line), except for 17 GeV flux, which was derived 
from only two photons. Because of its weak mag- 
netic field, less energetic synchrotron photons re- 
duce the absorption of the ICS component. The 
reduced absorption results in a small synchrotron 
flux, which further reduces the absorption outside 
the light cylinder, leading to a large, unabsorbed 
TeV fluxes. It also follows that the spectrum turns 
over at lower energy than the IG model prediction 
(dotted curve in the right panel; Harding 2001). 

It should be noted that the predicted ICS flux 
(above 100 GeV) represents a kind of upper lim- 
its, because it is obtained by assuming that all the 
magnetospheric synchrotron photons illuminate 
the equatorial region in which gap-accelerated 
positrons are migrating. Some of such VHE 
photons materialize as energetic secondary pairs, 
emitting the synchrotron component between a 
few keV and 100 MeV. Between 100 MeV and 
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30 GcV, the primary curvature component domi- 
nates, which represents the absolute flux (instead 
of upper limits). Some of such curvature photons 
have energies above 10 GeV and arc efficiently ab- 
sorbed to materialize as less energetic pairs, which 
emit synchrotron radiation below a few keV. Thus, 
the spectrum below a few keV also represents the 
absolute flux. 



4. Summary and Discussion 

To sum up, the self-consistent gap solutions ba- 
sically reproduce the observed power-law spectra 
below a few GeV for the three pulsars examined. 
The cut-off spectra between 10 GeV and 100 GeV 
strongly reflect the three-dimensional magnetic 
field configuration near the light cylinder; thus, a 
discrimination between IG and OG models (e.g., 
using GLAST) should be carefully carried out. 
If pulsations are detected above 100 GeV, it un- 
doubtedly indicates that the photons are emitted 
via ICS near the light cylinder, because VHE emis- 
sions cannot be expected in IG models. 

Since our analysis is limited within the two- 
dimensional plane formed by the magnetic field 
lines that thread the stellar surface on the plane 
containing both the rotational and magnetic axes, 
azimuthal structure is still unknown. Thus, we 
cannot present pulse profiles, phase-resolved spec- 
tra, or the polarization angle variations in this let- 
ter. Since the gap is most active in the outer part 
of the magnetosphere (unlike previous OG mod- 
els, which adopt the vacuum solution of the Pois- 
son equation and hence a uniform emissivity) , and 
since the photons will be emitted along the in- 
stantaneous particle motion measured by a static 
observer (unlike the treatment in the OG model 
of RY95, who assume a very strong aberration of 
light near the light cylinder) , it is possible that the 
predicted pulse profiles and so on are quite differ- 
ent from previous OG models. These topics will 
be discussed in the subsequent paper, which deals 
with the three-dimensional gap structure near the 
light cylinder. 

The pulsed TeV flux from PSR B1951+32 can 
be predicted to be above 5 x 10^° JyHz, pro- 
vided that a certain fraction (more than 30%) of 
the magnetospheric soft photons illuminate the 
equatorial region. However, if the poloidal field 



lines arc more or less straight near the light cylin- 
der, as demonstrated by Spitkovsky (2006, fig. 2) 
for an oblique rotator, the equatorial region may 
not be efficiently illuminated. In this case, the 
VHE flux will decrease from the current predic- 
tion. There is room for further investigation how 
to extend the present analysis into three dimen- 
sions, and to combine it with time-dependence 
three-dimensional force- free electrodynamics. 
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Advanced Research in Astrophysics (TIARA) op- 
erated under Academia Sinica and the National 
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Taiwan administered through grant number NSC 
95-2752-M-007-006-PAE. 
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Fig. 1. — Phase-averaged spectrum of the mag- 
netospheric emissions from the Crab pulsar for 
a; = 75°, /i = 4 X lO^o Gem^ and /i^ = 0.12. The 
thick dashed, dash-dotted, dotted and dash-triple- 
dotted lines represent fluxes into 93° < 9 < 97°, 
97° < 6* < 101°, 101° < 9 < 105°, 105° < 9 < 
109°; the thin dashed, dash-dotted, dotted ones 
109° < 9 < 113°, 113° < e < 117°, 117° < 9 < 
121°. The thick solid line represents the averaged 
flux in 75° < 61 < 103°. See Eikenberry et al. 
(1997) for IR-UV data; Knight (1982), Weisskopf 
et al. (2004), and Mineo et al. (2006) for X-rays; 
Nolan et al. (1993), Ulmer et al. (1995), Much et 
al. (1995), Fierro et al. (1998), Kuiper et al. (2001) 
for 10 MeV-20GeV; Borione et al. (1997), Tan- 
imori et al. (1998), Hillas et al. (1998), Lessard 
et al. (2000), and de Naurois et al. (2002) for the 
upper limits above 50 GeV. 
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Fig. 2. — Phase- averaged Vela spectrum for /i = 
4 X 10^0 Gcm^ ai = 75° and = 0.21. See 
Manchester et al. (1980) for optical data (open 
squares), Harding et al. (2002) for X-rays (filled 
triangles), Kanbach et al. (1994), Ramanamurthy 
et al. (1995b), Fierro et al. (1998), Thompson et 
al. (1999) for 7-rays (open circles). Left: Same 
figure as figure [TJ Right: Close up above 100 MeV 
to compare with IG (Harding & Daugherty 1993) 
and traditional OG (Romani 1996) models. The 
solid line is same as the left panel. 
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Fig. 3. — Phase- averaged spectrum of 
PSR B1951+32 for ^ 2 X lO^^Gcm^, a; 75° 
and /im = 0.39. Lines represent same quanti- 
ties as figure [2] unless notified. See Clifton et 
al. (1988) for IR upper limits (open squares), 
Becker & Triimper (1996), Safi-Harb et al. (1995) 
and Chang & Ho (1997) for X-ray data (open 
triangles), Ramanamurthy et al. (1995a), for 
100 MeV-15 GeV data (open circles), Srinivasan 
et al. (1997) for VHE upper limits. 



8 



